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We consider the femto-lensing due to a cosmic string. If a cosmic string with the 
deficit angle A ^ 100 [femto-arcsec] ^ 10~^^ [rad] exists around the hne of sight to 
a gamma-ray burst, we may observe characteristic interference patterns caused by 
gravitational lensing in the energy spectrum of the gamma-ray burst. This "femto- 
lensing" event was first proposed as a tool to probe small mass primordial black 
holes. In this paper, we propose use of the femto-lensing to probe cosmic strings with 
extremely small tension. Observability conditions and the event rate are discussed. 
Differences between the cases of a point mass and a cosmic string are presented. 
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I. INTRODUCTION 

Cosmic strings are line-like topological defects that are likely to emerge through phase 
transitions with spontaneous symmetry breaking Another possibility to produce cos- 

mic strings has been pointed out in the context of super-string theory and their properties 
are quite similar to those of field theoretic cosmic strings, except for the fact that inter- 
commuting probability between strings can be much lower than 1 (see [lo|-E3] for 
recent reviews). Usually, cosmic strings are characterized by the string tension /x. The ge- 
ometry around an infinite straight string is locally identical to that of fiat spacetime, but 
globally it corresponds to a conical spacetime with a deficit angle fl8| . 

A = ^. (1) 

where G is Newton's gravitational constant and c is the speed of light. One typical eflFect 
of the deficit angle is the formation of a double image of a light source located behind the 



string P^, Il9||. We call all such effects which are caused by gravitational fields of a cosmic 
string gravitational lensing due to a cosmic string. 

For ordinary field-theoretic strings, since the relevant energy scale is given by the sym- 
metry breaking scale Esb , the string tension can be estimated as Gji/c^ ^ E^^/Ep^ , where 
Ep\ is the Planck energy. For the Grand Unification scale, we have Esb ^ lO^^GeV or 
equivalently A ^ 10~^ . By contrast, the tension of cosmic strings from symmetry breaking 
along a supersymmetric fiat direction in the potential may have much smaller tension. In 
this scenario, there can exist two typical scales which determine the string tension. One is 
the supersymmetry breaking scale, which can be of the order of TeV, and the other is the 
cutoff scale ^ Ep\. Such strings may have the deficit angle given by 10~^^ < A < 10~^ ^ 



23|]. The effective four-dimensional tension of cosmic super-strings strongly depends on the 
details of the compactification and the infiationary scenario. In the case of the KKLMMT 
scenario jsl, the predicted tension of cosmic superstrings is expected to lie in the range 
10-12 <A< 10-6 0,0,0. 

So far, many attempts have been made to search for the signature of cosmic strings 
in various observations including the cosmic microwave background (CMB), gravitational 
waves and gravitational lensing. The CMB anisotropy spectrum has excluded cosmic strings 
with A > 10~^ from the dominant energy components of the universe Non-detection 



of g ravitational waves from cosmic string loops also rules out the string tension A ^ 10 



29|,l30| 



Gravitational lensing phenomena could in principle serve as more direct evidence for 



cosmic strings j31l-l34l| although none have been detected yet |35l-l37[| . A recent search for 
lensed galaxy pairs found no evidence for the presence of long straight cosmic strings with 
A > 7.5 X 10~^ out to redshifts greater than 0.6 [38]. Non-detection of characteristic 
variability of quasars due to the crossing of cosmic strings constrains lighter strings with 
10~^^ < A < 10"^ down to the level of = 0.01 j39|, where is the average density of 
cosmic strings in the units of the critical density. A more model-dependent but interesting 
bound on the local abundance of the strings with much lower tension 10~^^ < A < 10~^ 
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could be obtained by considering the variability of Galactic stars and pulsar timing j4( 
Forecasts for future constraints on the string tension include the B-mode polarisation 



of the CMB induced by straight strings with A < 10"^ [27|, l28|], weak lensing j4l|, ^ 



gravitational wave bursts, the stochastic background from string loops with A < 10 ^ j43| 



lensing at radio frequencies by loops with A ^ 10 [33|| and the 21 cm radiation from 
strings with A > 10-^0"^^ fl. 

So far, no way to probe cosmic strings with A < 10~^^ has been proposed. In this paper, 
we propose that femto-lensing events of gamma-ray bursts (GRBs) can be used to probe the 
cosmological abundance of cosmic strings with such small tensions. Though the expected 
image separation due to a cosmic string with A < 10~^^ is too small to be angularly resolved, 
the interference between the images could induce observable characteristic patterns in the 
energy spectrum of the lensed source objects |45|]. This effect, called femto-lensing |46|, 
was first proposed as a method to probe light compact objects like primordial black holes 
(PBHs) j46|-l49| and GRBs were considered as the target sources. The interference is expected 
when the time delay induced by lensing is comparable to the inverse of the gamma-ray 
frequency and femto-lensing of GRBs is sensitive to compact objects with a mass range 
10^7 g-10^° g. In the case of cosmic strings, as will be shown in the following sections, those 
with the deficit angle 10~^^~10~^'' can be probed by GRBs. 

Since the launch of the FERMI satellite, observational studies of GRBs have progressed 
significantly thanks to its unprecedented sensitivity. Recently, a constraint on the cosmo- 
logical density of compact objects in the mass range 5 x 10^'' g-10^° g was derived from the 
non-detection of femto-lensing events by the FERMI data in Ref. [50]. Likewise, we can 
expect a constraint on cosmic strings with tiny deficit angle 10~^^ < A < 10"^'', which can 
hardly be probed by other observations. 

The main purpose of this paper is to investigate the interference pattern in the energy 
spectrum of a GRB induced by a cosmic string and to determine how it differs from that for a 
compact object. The interference pattern for a compact object was discussed in Refs. (471. 1481 . 
Though typical gravitational lensing effects are often understood by using the geometrical 
optics approximation, the wave nature of the light is important in the case of femto-lensing. 
In the following sections, we discuss the observability of femto-lensing events due to cosmic 
strings by investigating wave propagation in a spacetime with a cosmic string [51]. 

This paper is organised as follows. In Sec. [TTl setting a lens system composed of the 
observer, a source and a straight cosmic string, we present the lensed waveform. A derivation 
of the lensed waveform using the Kirchhoff integral theorem is given in Appendix [Aj The 
observability conditions and the event rate are discussed in Sec. [TTTl In Sec. [IVl we point 
out some differences between the case of a cosmic string and that of a point mass. Sec. IVlis 
devoted to a summary. 



II. LENSED WAVEFORM 

In this section, we obtain the lensed waveform of a massless field in a cosmic string 
spacetime and see how the oscillatory behaviour in the energy spectrum can arise due to the 
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interference between the waves reaching the observer. 



A. Configuration of the Lens System 



We consider the spacetime with the deficit angle A, where A is defined so that the total 
angle around the cosmic string is given by 27r — A. Let us consider the configuration of 
the lens system specified by the following quantities: the distances from the observer to 
the string rfo, from the source to the string ds, from the observer to the source along the 
string direction dz and the angle ip which specifies the source position as shown in Fig. [H 
Hereafter, we assume (p = 0{A). 




O 



FIG. 1: Configuration of the lens system. The vertical straight line represents a cosmic string. The 
symbols O and S denote the positions of the observer and the source, respectively. 



An important quantity characterising the wave propagation in the lens system is the 
optical path difference I defined by 

_ dgdsA^ 
8D ' 

where 

D := ^{do + ds^ + dl. (3) 

The meaning of / is clarified by focusing on the (f = case, see Fig. [2l Then the distance 
from the observer to the source is given by (see the left of Fig. [2]) 



'dl + dl + dl + 2dods cos (^j^ =d(^1-^ + 0{A^)^ . (4) 

It is understood from Eq. that / approximately measures the shortening of the distance 
due to the presence of the string. 



dz 




o 



FIG. 2: Configuration of the lens system with ip = 0. 



For later use we also introduce an angle 



/3:= 



(5) 



2{do + ds)' 



Again the meaning of /3 becomes clear by setting ip = 0: then the half of the separation 
angle of the two images projected on the plane perpendicular to the string (see the right of 



Let us analyse the waveform after propagation from a stationary source to the observer. If 
we neglect the tiny tilt of the polarisation vector caused by the lensing, the local propagation 
equation for electromagnetic waves can be reduced to the wave equation for the massless 
scalar field 0: 



where uj is the angular frequency of the mode and (r, 6^, z) is a cylindrical coordinate system 
in which the cosmic string is on the z-axis. Since the spacetime with a deficit angle is locally 
flat, is just the flat Laplacian, but the range of the azimuthal angle 9 around the cosmic 
string is [0, 27r — A). 

Although we will not take cosmological expansion into account below, since the Maxwell 
equations in vacuum are conformally invariant, we can easily apply the results to cosmo- 
logical situations replacing Euclidean distances and the angular frequency uj with angular 
diameter distances and the redshifted angular frequency (1 + z)(jJ^ respectively, where z 
is the redshift at the intersection of the string and the line of sight and uj is the angular 
frequency at the observer. 



Fig. ED is /3 + 0(A2). 



B. Lensed Waveform 




(6) 
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Wave propagation in a locally flat spacetime with the deficit angle A has been fully 
studied in Ref. |5l|. Using appropriate approximations (see Ref. [si]]^ and Appendix El , we 
may solve the wave propagation to express the lensed waveform at the observer as 



(pido, ds, 4, ^; A, u) = F{w, y) 0o(c^o, 4, 4, ^; i^) , 



(7) 



where is the unlensed waveform (i.e. when A = 0) and F is the amplification factor, 
respectively, given by 



4>o{do, ds, 4, ^; ■= exp 



iuD 



F{w,y) := exp (1 + 2y) 



w 



dpds 2 



1 - -Erfc 
2 



(8) 



exp 



1 



-Erfc 
2 



:(! + ?/) 



(9) 



with A being an arbitrary amplitude, Erfc the complementary error function defined by 

2 



Erfc(x) := 



POO 

/ 

J X 



dt, 



and w and y the non-dimensional variables defined by 

2u!l ujdods/S? (pds 



w 



4cD 



y 



/3{do + 4) 



'a' 



(10) 



(11) 



It is worth noting that w gives roughly the ratio between the path difference and the wave- 
length, and y gives the source position Lpds on the source plane in the units of /3 {do + ds)- 
We are interested in the absolute square of F, which is observable as magnification or 
demagnification of the fiux. It is depicted as a function of w for several values of y in Fig. [31 

The oscillatory behaviour of |Fp in the frequency domain is of most importance for our 
purpose. In Fig. [3l the geometrical optics limit {w oc) of 



iFgeol' := lim \Ff = 



1 for \y\ > 1 

2 + 2co8{2wy) for \y\ < 1 



(12) 



is also depicted. Comparing |Fp with |FgeoP reveals that for most cases when y > 
the dominant mode in |Fp is the one that oscillates with period 6w = n/y. This is a 
consequence of the fact that the error function parts (those inside the curly brackets of the 
expression for F in Eq. (|9j)) only weakly depend on w. 

To summarise this section, we have analytically obtained the lensed waveform of the 
massless field in the cosmic string spacetime, Eq. ([7j). The numerical plots in Fig. [3] show 



^ We note that notations in this paper are totally different from Ref. [sil 
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that the absolute square of the amphfication factor |Fp osciUates in the frequency domain 
with the period Su given by 



Using the wavelength A, we obtain 

5uo A 2DX 



u Ayl dodsA'^y 



III. OBSERVABILITY CONDITIONS 



(14) 



In this section, we apply the above analysis to gamma-ray bursts from cosmological 
distances to determine the conditions under which the lensing effect due to a string can be 
observable. 



1. Detectable range of the deficit angle 

First, we evaluate the range of wavelength A within which a lensing event is detectable 
through the observation of photons with a fixed value of the deficit angle A. At least one 
whole phase of sinusoidal oscillation should be detectable in the spectrum. Hence, let us 
focus on one period of the sinusoidal oscillation around the angular frequency u) = cw/2l. 
Obviously, w should satisfy 2wy > n so that at least one period can be included in the 
spectrum. Thus, a necessary condition for an observable effect is given by^ 

A<8^/. (15) 

On the other hand, a lower bound for the wavelength is given by the energy resolution of the 
detector. The resolution of the detector must be sufficient to resolve the oscillation in the 
spectrum. We assume that the frequency resolution limit is given by auj with the typical 
value of a being (9(0.1). Then, we obtain the following inequality as another necessary 
condition: 

UJ Ayl 

Eventually, we find 

A<%/<a-iA « Aayl<\<Syl. (17) 

This implies that only a narrow range of wavelength, with the fractional width of order a, 
is usable in probing a given configuration of the lens system. We note that we are interested 
in the case (/; < A/2, or equivalent ly y ^ 1. 



^ As can be seen from Fig. [3l the gravitational lensing effect becomes only significant when w is larger than 
1, so we need ^ 1 ^ A < 47rL This inequality gives another upper bound for A. 
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If we take A = 10 ^ cm and 2do = 2ds = D = 10^^ cm with the observations of GRBs in 
mind, we have 

f Mods/D \ f ^ y ^ A f Mods/D \ f ^ y ...^ 

V 10^^ cm ; VlO-^V ^ 10-9 cm ^ ^ V lO^^cm ; \10-^y ' ^ ' 

Thus in the case of GRBs, the typical separation angle is evaluated as A ^ 100 femto arcsec, 
confirming the origin of the name "femto-lensing." Correspondingly, the detectable energy 
scale of the string tension is estimated as 



. 5 X 10»GeV (^y'V^^V", (19) 

VlO-^cmy V lO^^cm J ' ^ ^ 

where we have assumed I ^ X. 

In reality, the wavelength of GRB photons typically spans the finite range [Amin, Amax] = 
[10"^° cm, 10"'' cm], which roughly corresponds to the energy scale from IkeV to 1 MeV. 
Then the detectable range of the deficit angle is found to be 



2. Source radius 

So far, it has been assumed that a GRB can be treated as a point source. However, 
if the source radius is suflBciently large, the interference pattern will be smeared out. The 
upper limit for the source radius can be estimated by considering the y dependence of the 
amplification factor. For a fixed value of the square of the amplification factor oscillates 
with the period 5y = n/w diS di function of y. The corresponding length scale on the source 
plane to this period is given by 8y[5{do + rfs) = 7TAds/2w. If the source radius is larger than 
7rAds/2w^ the interference eflFect is smeared out and the oscillation in the spectrum cannot 
be observed. Since we are mainly interested in the situation / ^ A, we obtain w ^ An and 
find that the source radius must satisfy 

5. 6 X 10 cm — — — - — . (21) 

2w V47r/ V lO-isy V 1028cm y ^ ^ 



We estimate the source radius following Ref. [48[]. The appropriate linear source size 
is given by FcAt with V and At being the bulk Lorentz factor and the smallest variability 
time scale detected, respectively. Then we obtain 

rs - 3 X l{fcm ( ) [ \ • (22) 



100/ V 0.1ms 
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Once we fix the source radius Tg, tlie inequality f l2T]) represents tlie lower bound for the wave 
length A as follows: 

A > 1.5 X 10--cm ( -^^) ( ^1 ( . (23) 

Comparing this inequality with the range of wavelength ( fTSl) . we see that, although the 
finite source effect might be significant in some cases, the inequality ( I2T]) may be satisfied in 
many cases. We do not discuss the finite source effect in this paper. In a practical analysis, 
however, this effect may have to be taken into account (see Ref. |47| for the point mass case). 



3. Cosmic string motion 

In the previous discussions, we have assumed a static configuration of the lens system. 
However, if the relative velocity between the source-observer system and the cosmic string is 
too large, the string would pass through the region in which the lensing effect is significant 
before a detector collects enough photons. 

A straight cosmic string with relativistic vertical velocity v would pass through the region 
of interest in the time 

Mods/{do + ds) \ / lO^Qcm/s A ..4^ 
1028cm J \ V J ' 

A clear spectrum must be obtained within a time much shorter than this, otherwise the 
lens system cannot be regarded as static. For a large value of ^, only a limited number 
of bright burst events would be usable to detect femto-lensing events. This demand might 
be compared to those in the case of compact objects, which are believed to have a non- 
relativistic velocity dispersion. 



4. Cosmic string density 

Finally, we estimate how many strings in between the source and the observer are neces- 
sary in order that femto-lensing events can actually take place for a given number of GRBs. 
Here, we assume that the string network can be represented by a collection of straight string 
segments. 

Suppose a spherical volume of radius centered at the observer, contains A^cs straight 
strings of length L. Assuming the distances to the strings are ^ we obtain the angular 
scale of the length of a string as 2L/D. Therefore, the solid angle in which the lensing 
effect becomes significant is given by Ax2L/D. If there are no overlapping regions between 
the solid angles of neighbouring strings, the total solid angle given by all the strings is 
2AcsAL/L). Dividing this total solid angle by 47r, we obtain the event probability P for a 
single source as 
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Meanwhile, the average mass density is given in terms of A^cs by 

Introducing the density parameter of the straight strings Qcs = G Pcs/{3Hq)^ where Hq is 
the current Hubble parameter, we find 

p= — — . 27 

When we consider a source and a lens separated by cosmological distances, since D ^ cHq^^ 
we find P ^ Qcs- 

In reality, a GRB consists of many spike emissions, and cosmic strings are moving with 
the velocity v. Even if a femto-lensing event cannot be observed at the first spike emission, 
it might be observed at a subsequent spike emission during a single GRB event. This would 
in principle be possible if the GRB as a whole lasts longer than the crossing time scale given 
by (HM while there are spikes with timescales shorter than it. Still, those spikes may not 
individually contain enough photons to provide a clear spectrum; then we would need to 
collect a bunch of spikes. Let nb denote the mean number of such qualifying bunches within 
one GRB. Then, the event probability P for a single GRB has to be multiplied by the factor 
riiy. Assuming f^cs ^ 0.001 and nb ^ 10, we can expect roughly one femto-lensing event 
among 100 available gamma-ray burst events. This roughly corresponds to one detection 
per year with FERMI satellite if most of the observed gamma-ray bursts can be used for 
the femto-lensing search. Even if we do not observe such an event, we may obtain an 
observational limit on f^cs- Further detailed analysis of GRB spectra is needed for a more 
precise estimation of the event rate and giving observational limits on f^cs- 



IV. DISTINCTION BETWEEN COSMIC STRING AND POINT MASS LENS 



Cosmic strings are not the only candidates for a femto-lensing object, as point masses 
would also cause a similar eflFect. It is necessary to study the differences between these two 
cases in order to determine how they can be distinguished observationally. Below, we shall 
review the femto-lensing caused by a point mass and highlight the diflFerence between this 
and the case of a string. 

Let us consider the point mass lens system with distances from the observer to a lens 
L^L, from the lens to a source L^ls and the observer to the source Ds = Di^ + Di^s- Let 
T^p denote the distance between the source position and the intersection of the source plane 
and the line which connects the observer and the lens. For the point mass case, we have the 
following expression for the amplification factor [52|, |5J] : 



iw 



2 ' ' 2 J 



(28) 
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where F and iFi are the gamma function and the confluent hyper-geometric function, re- 
spectively, and 



AojGM 



4GMDi^sDs 



(29) 



Considering the rjp = case, we can find the following expressions for the path difference 
Ip and half of the separation angle /?p: 



2GM 

r.2 ' 



4G'MDls 



(30) 



These /p and j3p correspond to / and j3 in the cosmic string case. Using these quantities, we 
can express Up and Wp as 



2uj L 



Vp 



Vp 



/3p(Dl + Dls)' 



(31) 



One can see that these expressions for and are similar to those for y and w in the 
cosmic string case f fTTI ). Therefore, we identify y^ and as y and respectively, and omit 
the subscript p hereafter. 

In the geometrical optics approximation {w oo), we obtain 



geol 



Im+I + IM-I + 2-^|/x+//_| cos(0_ - ^+), 



where 



±- 



1 



+ 



±2 



24 
1 

2Z 



In 

In \x- 



71 

2'' 



(32) 

(33) 

(34) 
(35) 
(36) 

(37) 

cos(^_ - ^+) = ^m{wT{y)). (38) 

It is meaningful to compare r(|/) with 2y because we have Eq. flT^ for the cosmic string 
case in the geometrical optics approximation. As is shown in Fig. HI riy) is close to 2y for 

Before comparing the cosmic string case with the point mass case, we compare the am- 
plification factor with that of the geometrical optics approximation for each case. As is 



^_ — Qj^ can be rewritten as 
Therefore we obtain 



-2/ V?/^ + 4 + In 

2 + 4 _ ^ 



2 ^^^^^ - 2 
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FIG. 4: T{y) and 2y. 



Point mass, y-ljl 



In An 6n Stt IOtt 12;: Un I6n ISn 20n 




2n 4n 6n Stt IOtt 12;: Un I6n IStt 20;: 



w 

FIG. 5: and |Fg eo| cis functions of w for the point mass case with y — 1/2. 



shown in Fig. El the spectrum quickly gets closer to the geometrical optics approximation 
in the point mass case. Though we show only the result for y = 1/2, this behaviour is 
a common feature for any value of y. On the other hand, this is not the case when we 
consider the cosmic string, as is shown in Fig. [6l While the phase of the oscillation in the 
spectrum seems to be well approximated by the geometrical optics approximation, there are 
remarkable deviations in the amplitude. In this sense, we can conclude that the wave effect 



14 



Cosmic string , y=ll2 



In 4n 6n Stt lOn I2n Un I6n ISn 20n 




In 4n 6n Stt lOn I2n Un I6n ISn 20n 



FIG. 6: |Fp and |i^geoP as functions of w for the cosmic string case with y — 1/2. 



is more important for the cosmic string case. 

Finahy, let us look at the differences between the cosmic string case and the point mass 
case. In order to highlight the differences, we focus on the following quantity: 





F\ 


2 


( ^geo 





where (|FgeoP) is the 1 cycle average of the square of the amplification factor in the geomet- 
rical optics approximation. Using this procedure, we can extract the oscillating part of the 
spectrum. This quantity is plotted as a function of w for each case in Fig. [71 

We find two remarkable differences between the two cases. One is the 7r/2 phase difference, 
the origin of which is the — 7r/2 shift in Eq. f l35l ). This shift is caused by the caustic of one 
of the ray bundles in the case of the point mass lens. This does not happen in the 
cosmic string case. The other difference is that the oscillation amplitude of |Fp/(|FgeoP) 
can exceed 1 in the cosmic string case, but cannot in the point mass case. The reason 
why it does not occur in the point mass case is obvious from Eq. f l32l ) and the goodness of 
the geometrical optics approximation. These differences might be used to identify the lens 
object if we succeed to detect the oscillatory behaviour in the spectrum of a GRB. 



V. SUMMARY 

Femto-lensing events of gamma-ray bursts (GRBs) due to cosmic strings have been 
studied. The detectable range of the deficit angle A was derived taking the wave na- 
ture of gamma-rays and the energy resolution of a detector into account. Assuming the 
range of observable wavelengths to be 10~^°cm-10~''cm, we obtain the detectable range as 



15 



Cosmic string(CS) and point mass(PM), y=ll2 
In An 67r Stt IOtt lln lAn IGtt 




w 



FIG. 7: |Fp/(|FgeoP) — 1 as a function of w for the cosmic string case and the point mass case 
with y — 1/2. 

20-19 < /\ < 10"^''. Observabihty conditions associated with the source radius and the 
relative motion of the lensing system have also been discussed. The relative motion of a 
cosmic string gives a limitation on the duration of photon counting. This limitation may be 
tighter than that in the case of a compact lens object because of the relativist ic motion of a 
cosmic string. A limited number of bright burst events can be candidates for a femto- lensing 
event if the relative velocity of a cosmic string is comparable to the speed of light. The event 
probability P for a single GRB event is roughly estimated as P ^ ^b^cs, where rib and f^cs 
are the mean number of available bunches of spike emissions in a GRB and the average 
density of cosmic strings in units of the critical density, respectively. 

Two typical differences between the lensed spectrum in the case of a point mass lens and 
that of a cosmic string have been pointed out. One of these is a phase shift in the spectrum 
oscillation. In the point mass case, one of two ray trajectories experiences a caustic which 
causes a — 7r/2 phase shift, and this shift appears in the lensed spectrum. The same does 
not occur in the cosmic string case. The other difference is in the oscillation amplitude of 
the spectrum. In the point mass lens case, the oscillation amplitude does not exceed the 
mean value, while it can exceed the mean value in the cosmic string case because of the wave 
effect. These differences might be used to identify the lens object if we succeed to detect 
the oscillating behaviour in the spectrum of a GRB. 
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Appendix A: Derivation of the Lensed Waveform 



We give an overview of the derivation of Eq. ([71) using the Kirchhoff integral theorem [53 
This approach for the point mass lens case can be seen in Ref. [5J]. A diflFerent approach is 
used in Ref. [5l(] to derive Eq. ([71). We use the unit c = 1 in this Appendix for notational 
simplicity. First, let us introduce the coordinate system ^ = ((^, C z) where the observer and 
string are located at (0, —do^O) and C = C = 0, respectively. For convenience, we remove 
the region of the deficit angle so that the source is just split in two as shown in the left of 
Fig. [HI In this coordinate system, the source position can be expressed in two ways. Let 




identical 




FIG. 8: Relation between the coordinate system z) and the configuration of the lensing system 
is shown in the left figure. Another identical description which has been used in the text is shown 
in the right figure. 



i^sXs^dz) denote the two expressions for the source position, where > 4 • Then, we 
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find 



4 = ids sin ( — ± ^ ) , 



Cs = dscos[—±(p] . 



(Al) 
(A2) 



Our purpose is to obtain an approximate solution for the wave equation given by 

^^^+uj^)^ = -4nAS{^-^,). (A3) 

The strategy is summarised as follows. First, we assume that the geometrical optics ap- 
proximation is valid in the domain ( > 0. Then, we calculate the waveform at the observer 
by applying the Kirchhoff integral theorem to the domain ^ < 0, where the waveform of 
the geometrical optics approximation on the C = plane is used as a boundary condition. 
To evaluate the Kirchhoff integral, we use appropriate approximations associated with the 
small quantities ip ^ A ^ e and 1/{ujD). 

Using the geometrical optics approximation, the waveform around the C = plane is 
given by 



exp{iujDf) for ± ^ > 0, 



(A4) 



where 



= - ^ty + (c - ctr + - d.y. (A5) 

Applying the Kirchhoff integral theorem to the domain C < and neglecting the contribution 
from the boundary at infinity, we obtain the following integral for the waveform at the 
observer: 



bo := 0(0, -do, 0) = 



1 

47r 



dz 



d_ 



~d7 



d 



(A6) 



C=o 



where 



Do 



'e + {C + doY + z^- (A7) 

In order to obtain an approximate form for this integral, we keep terms up to the order 
of in the phase part and up to the leading order in the amplitude. It is intuitively obvious 
that only the region ^ < eD can give significant contribution to the integral. This fact 
can be justified by applying the stationary phase approximation to the integral over ^ with 
1/{(jjD) <C 1. Taking the above discussion into account, the integral can be approximated 
as 



-iujA 
47r 



dz 



o 



X 



exp 



+ exp 



-lU) 



-lUJ 



D3DI 

2 



+ 



exp [iu {D3 + D^)] 



2D3 {Ds + D5) 

_2 



2D3 {D3 + D5) 




d^exp 



d^exp 



iujDs + D5) 



D3+D5 



,(A8) 
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where 



= ^dl + {z- 4)2, 
D, = 



Defining a and by 



a = 



^5^ 



(A9) 
(AlO) 

(All) 
(A12) 

(A13) 

The integral with respect to z can be evaluated by using the stationary phase approx- 
imation. Since the dominant contribution in the phase part comes from i(jj{D^ + D5), we 
write the Eq. (|A8j) in the following form 



we can evaluate the integral of ^ as 

d^exp [ta{±^ - b^f] = (1 - ^Erfc (±V^b^^'\ . 



9o 



dzX{z) exp[iw(L>3 + D^)], 



(A14) 



where 



X{z) 



-iuA / (io _|_ c?s 



X 



47r \DsDl DID, 



exp 



-lU) 



+ exp 



2Ds {Ds + D^) 

_2 



-lU) 



2D3 {Ds + D5) 



y^(l-iErfc(»)) 

7 



(A15) 



Then, we consider the stationary phase approximation with the phase iuj{D^ + D5). Since 
the stationary point is given by 



az do + rfs 



(A16) 



Eq. f lA14p can be approximated by 



/oo POO r 

dzX{z) exp[iuj{Ds + D^)] ^ X{z') / dzexp iu <D + 
-00 J —00 L L 



(do + dsY 



{z - z'f 



2dodsD^ 



Xiz'), 



2widodsD^ 



exp[iuD]. 



coido + dsY 

Calculating X{z') and simphfying the expression, we finally get Eq. 



(A17) 
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